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Abstract. In recent years, there has been increasing interest in consecutive pattern avoidance 
in permutations. In this paper, we introduce two approaches to counting permutations that 
avoid a set of prescribed patterns consecutively. These algoritms have been implemented in 
the accompanying Maple package CAV, which can be downloaded from the author's website. As 
a byproduct of the first algorithm, we have a theorem giving a sufficient condition for when 
two pattern sets are strongly (consecutively) Wilf-Equivalent. For the implementation of the 
second algorithm, we define the cluster tail generating function and show that it always satisfies 
a certain functional equation. We also explain how the CAV package can be used to approximate 
asymptotic constants for single pattern avoidance. 



1. Introduction 

Let (7 = (Ti • • • (Tfc be a sequence of k distinct positive integers. We define tfie reduction p{a) 
to be the length k permutation we get by relabehng the elements of a with {1, . . . , fc} so that 
they retain the same order relationships they had in a. For example, p(5386) — 2143. For a 
permutation p, we will also write \p\ for the number of elements in the permutation. Let m and 
n be positive integers with m < n, and let p G Sm and cr = cti • • • (t„ ^ Sn- We will say that a 
contains p consecutively if p{ai ■ ■ ■ cFi^rn-i) = P for some i where l<i<n — m + l. Otherwise, 
we say that a avoids p consecutively. Similarly, if _B is a set of permutations, then we say that a 
avoids B consecutively if for every p €z B, a avoids the pattern p consecutively. For example, the 
permutation 123654 £ Sq contains the permutation pattern 1243, since p(2365) = 1243. However, 
the permutation 12453 £ 55 avoids the pattern 1243 consecutively. 

In general, we are interested in counting permutations that avoid a pattern (or a set of patterns). 
Given a set of patterns B, let asin) be the number of length n permutations that avoid B 
consecutively. If B consists of only a single pattern p, we may write ap{n) instead, and if no 
ambiguity would arise, we may just write a{n). For a given set of patterns S, we would like to 
find the exponential generating function 

°° z" 

(1) = 

If no ambiguity would arise, this may also be denoted by A{z). In addition, we define a more 
general exponential generating function 



(2) PBiz,t) = 



k,n>0 



where bk^n is the number of length n permutations that contain exactly k occurrences of the 
patterns in B. Again, we may write P(z,t) if the set B is clear. We will also define a{n,t) = 
Y.k>Q^k,nt''. Note that P(z,0) = A{z) and a(n, 0) = a{n). 
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In addition, we will say that two sets of patterns B and B' are consecutively Wilf-Equivalent 
(sometimes written c- Wilf-Equivalent) if Ab{z) = Ab'{z). We will also say that B and B' are 
strongly c- Wilf-Equivalent if PB{z,t) = PB'{z,t). Since this paper deals solely with consecutive 
patterns, the word "consecutive" will be omitted in most instances. For the rest of this paper, 
the reader should assume that all mentions of containment, avoidance, and Wilf-Equivalence are 
consecutive. 



In recent years, there has been an increasing amount of research done on consecutive pattern 
avoidance in permutations. One of the early papers by Elizalde and Noy ( [S]) finds generating 
functions A{z) and P[z,t) for certain cases of single pattern avoidance. Using various techniques, 
additional generating functions for specific single patterns and multi-pattern sets have been found 
in [l]|3l|8l[3. In particular, our approach will resemble the cluster method approach in [3]- 

The results in this paper utilize an extension of the Goulden- Jackson cluster method ([SI [TU]V 
We restate some of the terminology and notation here. 

Let i? be a set of patterns. Without loss of generality, assume that B contains no trivial re- 
dundancies (i.e., there are no pi,P2 G B with pi ^ p2 such that pi contains ^2)- We say that an 
ordered pair (p; . . . , [im, jm]]) is a length k cluster if it satisfies the following: 



(a) p £Sk 

(b) Zl = 1, 

jm — k, and in < in+i < jn for 1 < ?i < m — 1 (i.e., each interval overlaps with 
the neighboring interval, and the intervals cover p) 

(c) p(pi„ • • -Pjn) € B for all 1 < rt < TO. 



Let Cfc be the set of clusters of length fc, and for a cluster w — {p; . . . , [jm,im]]), define 

weight{w) = (t — 1)™, where t will be the variable used to track occurrences. Let C{k) — 
^ weight[w) . From an adaptation of |Q to the present context of an "infinite" alphabet and 

exponential generating functions, we have: 

Theorem 1. 

(3) P{z,t) 



fe>i 



This theorem will essentially be the basis for our algorithm and our results. 

So far, generating functions have been found for specific single patterns and multi-pattern sets 
and for certain single pattern families where some specific structure can be exploited. In this 
paper, we will outline two algorithms to calculate a(ri) more efficiently, and both algorithms have 
been implemented in the accompanying Maple package CAV. The Maple package can be down- 
loaded from the author's website. As a result of the first algorithm in Section [21 we get a theorem 
for proving when two pattern sets are strong c- Wilf-Equivalent. During preparation of this paper, 
the author learned that this result was also independently proven by Khoroshkin and Shapiro in 
[3 by slightly different means. To establish the much faster second algorithm in Section [31 we 
define a new generating function which we refer to as the cluster tail generating function. We show 
that this generating function always satisfies a certain functional equation and give a constructive 
approach to finding it. This functional equation is then used to compute values for a(n) much 
more quickly. We use our algorithm to give some asymptotic approximations in Section [4l We 
conclude with Section [5l by sharing some new conjectures we have based off of experimentation 
with our CAV package. Beyond the theorems and results in this paper though, we hope that the 
CAV Maple package will be a useful tool for others in studying consecutive pattern avoidance in 
permutations. 



COMPUTATIONAL APPROACHES TO CONSECUTIVE PATTERN AVOIDANCE IN PERMUTATIONS 3 



2. Consecutive Pattern Avoidance via Clusters 

Let i? be a set of fixed patterns that we would like to avoid (consecutively). We may assume 
that B contains no redundancies (i.e., there does not exist pi,P2 G B with pi ^ p2 such that pi 
contains p2)- Again, a{n) will be the number of length n permutations avoiding B. 

From the Goulden- Jackson cluster method (Theorem [ij , we can get the equation 

(4) P(z, t)^l + zP{z, t) + P(z, t) J2 C{k) ^ 

fc>i 

and by extracting the coefficients of z", we get the following recurrence: 

(5) a{n,t) = na(n - 

k=l ^ ^ 

Additionally, consider a fixed p £ B and let m — \p\. Let Ck[p] — {(tt; [[ii, ji], . • . , [«r, jV]]) G Ck '■ 
p{TTi^ • • • TTj^) — p}, the set of length k clusters ending in the pattern p. Let Ck 11 be 

the clusters in Ck[p] with the last m terms {xi, . . . ,0;™}, where xi < X2 < ■ ■ ■ < Xm- Similarly, 
define 



(6) C{k,p) — weight{w) 

(7) C(fc,p; [xi, . . . = ^ weight{w). 

iueCfc[p;[a:i,...,a:,„]] 

If B contains only one pattern, these may be denoted by C(fc) and C(fc; [xi, . . . , Xm\), respectively. 
We will use Equation ([5]) to count permutations avoiding B. 



2.1. General Algorithm. Computationally, the difficulty in using Equation ([5]) lies in calculat- 
ing C{k) quickly. One way to do this is to create a recurrence for C'(k,p; [xi, . . . ,x\p\]) for each 
p€B. 

We can do this as follows: for a given cluster w, let pi and p2 be the last marked pattern and the 
second to last marked pattern in w, respectively. Let j be the length of the overlap of pi and p2 in w, 
i.e., the tail of length j of p2 coincides with the head of length j of pi. We want to "chopoff" the last 
\pi\ — j terms of w and apply the reduction to get a shorter cluster, say w' e Ck'[p2] [xi, ■ ■ ■ ,a^|p2|]L 
which ends in the pattern p2- Then, weight{w) = weight{pi) ■ weight(w'). 

Additionally, once we "chop off" the tail of pi and apply the reduction to get a shorter cluster w', 
we actually know what the last j terms of w' = w'l - ■ ■ w'j^, will be. For each term w[ with — + 
1 < i < the reduction forces w[ to be ■Wi — {^ of terms in w "chopped off" that were less than Wi). 
Thus, to compute C'{k,pi; [xi, . . . , a:|pj^|]), we need to sum over all possible ways to "fill out" the 
rest of the terms in the final p2 pattern of w' . We also need to sum over all possible choices of 
P2 & B and all possible ways that the tails of this p2 overlap with the heads of the final pi pattern. 

In summary, the number of length n permutations avoiding set B can be found by, first, 
generating a cluster recurrence for C'{k,p; [xi, . . . ,x\p\]) for each p € B. Next, use the recurrence 

n 

a{n) = na{n — 1) + X) (fc)^(^)'^("- ~ ^) using the base cases a{0) — a{l) = 1 and a{n) = if 

k=l 

n < 0. Use the cluster recurrences to compute C{k) as needed: 
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(8) C{k) = J2C{k,p) 

peB 

(9) C{k,p) = Yl C{k,p;[xi,...,xm]) 

l<Xi<X2<---<Xm^k 

Also recall that if w = (tt; . . . , [im,jm]), then weight{w) = [t — 1)™ will keep track of 

occurrences of patterns with variable while setting t = Q and using weight{'w) = (—1)™ would 
count only the permutations that avoid the designated pattern set B. 

2.2. Example. Let B = {2143}. Let w = (tt; [[ii, ji], . . . , [im, jm]]) be a length k cluster and 
{a;i, . . . , X4} be the last 4 terms of w with xi < . . . < X4 (i.e., 7rfe_3 = X2, 'Kk-2 = xi, nk-i = X4, 
and TTfe = X3). Then, the second to last pattern must also be a 2143 pattern and can have an 
overlap of length 1 or 2 with the last pattern. 

If the overlap is of length 2, let vr' = p{7:i ■ • ■TTk-2) and w' = (tt'; [[ii, ji], • . • , [jm-i, im-i]]), 
the cluster found by "chopping off" the tail of the final bad pattern in w and then canonically 
reducing. Now let {yi, ■ ■ ■ ,1/4} be the last 4 terms of w' with yi < ... < 2/4 (i.e., 7r^_3 = y^, 
7r^_2 = ^fe-i = 2/4: and Tr'f. = ya). Notice that the terms "chopped off" from w were X4 and X3. 
Since both of these are larger than both X2 and Xi , applying the reduction does not change their 
values. Thus, 2/4 = X2 and y3 = xi. Summing over all possible tails for w' and accounting for the 
last pattern that was removed from w, we get 

(10) J2 weight{2U3)-C{k-2;[yi,y2,y3,y4])- 

l<yi<...<yi<k-2 

y3=xi 

If the overlap is of length 1, let tt' = p(7ri • • • 77^-3) and w' = (tt'; [[ii, ji], . . . , jm_i]]), 
since the tail that gets "chopped off" has 3 terms. Again, let {2/1, . . . , 2/4} be the last 4 terms of w' 
with 2/1 < • • • < 2/4 (i-6-, "^'k-s, = 2/2, ^fe_2 ^ 2/1' '^fe-i ^ 2/4' ^^"^ ''^'k ~ V^)- '^be terms "chopped off" 
from w are xi, X4, and X3. Since exactly one term less than X2 (only xi) was removed, applying 
the reduction would reduce X2 by 1. Thus, 2/3 = 0:2 — 1. Summing over all possible tails for w' 
and accounting for the last pattern that was removed from w, we get 

(11) iveight{2U3)-C{k-3;[yi,y2,y3,y4])- 

l<yi<...<Vi<k-3 
2/3=2:2-1 

We combine the two possibilities along with the base cases to get the recurrence. 
For /c < 4: 

C{k; [a;i,a;2,a;3,a;4]) = 

For k = A: 

C{k;[xi,X2,X3,X4]) = weight{21'l3) 

For A: > 4: 

(12) C{k;[xi,X2,X3,X4]) = ^ u;ei5f/it(2143) ■ C(fc - 3; [2/1,2/2,2/3,2/4]) 

l<J/l<...<l/4<fc-3 

+ YI weight{2U3) ■ C{k - 2; [2/1, 2/2,2/3, Vi]) 

l<j/i<...<I/4<fe-2 

V4=X2 
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Using this recurrence, we can compute C{k) for any value of k and compute a{n) using Equation 
([5]). To keep track of all occurrence of 2143 with the variable t, let weight{214:3) = t — 1. To only 
count permutations that avoid 2143, set i = so that weight {2143) = —1 for the above recurrence. 

2.3. Results for c-Wilf-Equivalence. Even though Section lTll is algorithmic in nature, it yields 
a strong theoretical byproduct. The cluster recurrence generated by the pattern set B totally de- 
termines a(n). In fact, it also totally determines P{z,t). However, the "overlapping" relations 
between the patterns in B totally determines the cluster recurrence. 

More specifically, let B be the set of patterns we want to avoid, and let tt, cr G _B where m — \tt\ 
and n — \a\. Note that n and cr are not necessarily distinct. Suppose that p((T„_j+i ■ ■ - an) ~ 
p{'!Ti • • • TTj) (the tail of a and the head of tt has an overlap of length j). Then, define the following 
sets: 

(13) OverlapMap{a,Tr,j) = {(tti, a-„_j+i), (7r2, (t„-j+2), ■ • ■ , (tt^ , cr„)} 

(14) OverlapMaps{a, tt) = {OverlapMap{a, tt, j) : p{an-j+i ■ ■ ■ cr„) = p{Tri ■ ■ ■ Wj)} 



For example, in Section 12.21 the pattern 2143 has self-overlaps of length 1 and 2. For a 
length 1 overlap, we have 0-yer/apMap(2143, 2143, 1) — {(2,3)}. This combined with the length 
of the pattern, which is 4, and the length of the original cluster, denoted by k, completely 
determines the first summation in Equation (jl2p . Similarly, for a length 2 overlap, we have 
OverlapMap{2143, 2143, 2) = {(2, 4), (1, 3)}. Combining this with the length of the pattern, again 
4, and the length of the original cluster, again fc, completely determines the second summation 
in Equation ([H]). Thus, Ower/apMaps(2143, 2143) = {{(2, 3)}, {(2, 4), (1, 3)}} and |2143| = 4 
completely determines the cluster recurrence. 

Therefore we have the following result based off of our algorithm: 

Theorem 2. Let B and B' be two sets of patterns with \B\ = \B'\. Suppose there is some labeling 
of the elements (patterns) in sets B and B' , say B — {pi, . . . ,pk} and B' = {p'l, . . . such that 

\pi\ = \p'^\ for 1 < i < k, and OverlapMaps{pi,pj) — OverlapMaps{p[,p'j) for all 1 < i,j < k. 
Then, B and B' are strongly c- Wilf-Equivalent. 

Proof. The cluster recurrence was uniquely determined by how the patterns overlapped (which 
terms from one pattern overlapped with which terms of another pattern) and by how they reduced 
after "chopping" the last pattern from the current cluster. The possible ways that two patterns 
can overlap are encoded by OverlapMaps and the effect of the reduction p is determined by how 
the patterns overlapped and the length of those patterns. □ 

This result was also independently discovered by Khoroshkin and Shapiro, [7]. 

Using this result, it is possible to classify c-Wilf-Equivalences in some cases. For example, it 
is possible to classify single pattern avoidance for single patterns of length 3, 4, and 5 since all 
the potential equivalences that occur can be demonstrated using Theorem [5] Using the same 
approach, we can nearly classify all single patterns of length 6. All that remains are four possible 
strong c-Wilf-Equivalences that appear true but cannot be rigorously proven through our means. 
They are the following: 

(1) The pattern 123546 appears to belong to the strong c-Wilf-Equivalence class {124536, 125436}. 

(2) The pattern 123645 appears to belong to the strong c-Wilf-Equivalence class {124635, 126435}. 
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(3) The patterns 132465 and 142365 appear to be strong c-Wilf-Equivalent. 

(4) The patterns 154263 and 165243 appear to be strong c-Wilf-Equivalent. 
The four cases have been experimentally verified for up to length 18 permutations. 

2.4. Maple Implementation. The algorithm from Section 12.11 has been implemented in the 
Maple package CAV. Using that algorithm, you can find the sequence a(l), . . . , a{n) avoiding a set 
of patterns B by calling the procedure CAV(B,n) , where the patterns in B are represented as lists. 
For example, for n = 10 and S = {123,321}, trying CAV({[1,2,3] , [3, 2, 1] } , 10) ; returns the 
output: 

[1, 2, 4, 10, 32, 122, 544, 2770, 15872, 101042] 

To keep track of the occurrences of patterns from B, use the procedure CAVt (B , n , t ) . For example, 
trying CAVt ({ [1 ,2 ,3] , [3, 2 , 1] } ,6 ,t) ; returns the output: 

[1, 2,4 + 2t, 10 + 12t + 2t^ 58t + 28<^ + 32 + 2^^ 300t + 236^^ + 122 + 60t^ + 2^^] 

Also, most of the main procedures in the Maple CAV package have an optional verbose setting. 
For example, for the verbose outputs, try CAV({[1,2,3] , [3,2,1] }, 10, true) ;. 

To generate the cluster recurrence only (encoded in a data structure that we call a scheme), use 
the procedure SCHEME(k,B,x,y,t). For example, try SCHEME(k, { [1 , 2 , 3] , [3,2 , 1] } ,x,y ,t) ; . 
The overlap maps between two patterns can also be found using OverlapMaps (pi ,p2) , where you 
are checking for overlaps between tails of pi with heads of p2. For example, try OverlapMaps ([2,1,4,3] , 
[2,1,4,3]) . 

To (attempt to) classify pattern sets of m patterns with each pattern length n, we can compute 
a{N) (for some fixed value N) for each of these pattern sets, and if the a{N) values coincide, 
try to apply Theorem [2] This has been implemented in the procedure Wilf Eqm(n,N,m) . For 
example, try Wilf Eqm(5 , 12 , 1) (or for the verbose output, Wilf Eqm(5 , 12 , 1 , true) ) to (rigor- 
ously) classify c-Wilf-Equivalence for all single patterns of length 5. An additional byproduct 
of Theorem [5] is that all instances of c-Wilf-Equivalence in single length 5 patterns are actually 
strong c-Wilf-Equivalence. The 25 c-Wilf-Equivalence classes can be found on the paper's website. 

Similarly, we can use the Wilf Eqm procedure to discover the following: 

Proposition 1. Let Bi and B2 both he sets containing two patterns of length 3. Then Bi is c-Wilf- 
Equivalent to B2 if and only if they are trivially equivalent by reversal and/ or complementation. 

Proof. Run "WilfEqm(3, 10,2,true) ;" using the CAV Maple package. □ 

Proposition 2. Let Bi and B2 both be sets containing two patterns of length 4. Then Bi is c-Wilf- 
Equivalent to B2 if and only if they are trivially equivalent by reversal and/ or complementation. 

Proof. Run "WilfEqm(4, 10,2,true) ;" using the CAV Maple package. □ 

Proposition 3. Let Bi and B2 both he sets containing three patterns of length 3. Then Bi is 
c-Wilf-Equivalent to B2 if and only if they are trivially equivalent by reversal and/or complemen- 
tation. 

Proof Run "WilfEqin(3, 10,3,true) ;" using the CAV Maple package. □ 

Similarly, nearly all c-Wilf-Equivalences could be classified for sets containing three patterns of 
length 4. Four pairs of sets appear c-Wilf-Equivalent but cannot be proven through our means. 
They are the following: 
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(1) The pattern sets {1234, 1243, 1342} and {1234, 1243, 1432} appear to be strongly c-Wilf- 
Equivalent. 

(2) The pattern sets {1234, 1243, 2341} and {1234, 1243, 2431} appear to be strongly c-Wilf- 
Equivalent. 

(3) The pattern sets {1324, 1342, 1423} and {1324, 1423, 1432} appear to be strongly c-Wilf- 
Equivalent. 

(4) The pattern sets {1324, 1423,2341} and {1324, 1423,2431} appear to be strongly c-Wilf- 
Equivalent. 

The four cases have been experimentally verified for up to length 14 permutations, and the rest 
of the classification can be found on the paper's website. 

3. Consecutive Pattern Avoidance via the Cluster Tail Generating Function 

Computationally, the cluster recurrence is faster than the naive approach of checking every 
single permutation, but the approach is still very inefficient. For a fixed length fc, not every com- 
bination of tails gives rise to a possible cluster. For example, ii B = {123}, the only possible 
underlying permutation in a length 9 cluster is 123456789. The only possible tail is 789, but using 
the recurrence, we essentially try all (g) possible tails. Each such possible tail gives its contribution 
of only after it has recursed down to the base cases of fc < 3. 

We can, however, gain a substantial speed-up by considering a more complicated generating 
function. For a fixed pattern p ^ B with length m, the cluster tail generating function will be 
defined as: 

(15) F{k,p;[zi,...,z^]) ^ ^ C(fc,p; [xi, . . . • • -z^'" 

If B is a single pattern set, this may also be denoted as F{k; [zi, . . . , Zm])- Otherwise, we also 
define: 

(16) F{k;[zi,...,z^]) = ^F{k,p;[zi,...,Zm]) 

Note that F{k,p; [1, . . . , 1]) = C{k,p) and F{k; [1, . . . , 1]) = C{k). In fact, we can always find 
a functional equation for F(k,p; [zi, . . . , Zm]) of a certain form. We can then combine this with 
Equation ([5]) to more quickly compute a{n). We begin with an illustrative example and then 
present the general algorithm. 

3.1. Example. Let B = {132} and suppose we want to only count permutations that completely 
avoid 132. We will set t = which gives us weight{132) = —1. We then can find a functional 
equation for F{k; [zi, 22,^3]) as follows. Using the procedure SCHEME in the Maple package CAV, 
we can get the following cluster recurrence: 

C{k;[xi,X2,X3]) = - ^ C(/c - 2; [yi, 2/2,2/3]) 

i<yi<y2<y3<k-2 
y2=xi 

= ^ C'(fc-2;[yi,a;i,y3]) 

l<yi<xi 
xi<y3<k~2 
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with the base cases C(fc; [xi,X2, x^]) = if fc < 3 and C(fc; [xi, X2, x^]) = — 1 if fc = 3. Substituting 
into Equation (jl5p and applying the finite geometric series formula as needed, we get: 



^2 ^3 



^2 -'S 



i^(fc; [zi,Z2,Z3]) = ^ C(k\\xi,X2,x-:^zl'' 

l<a;i<a^2<a;3</c 

k-2 k-l k 

= -E E E E C{k-2;[y,,x,,y,])zr 

a;i = l 2:2=2:1 + 1 a;3=a;2 + l l<yi<xi 
xi<y3<k-2 

fc-2 k-l k 

= -E E E C{k-2;[y,,x^,y^M^zr ^ 

Xx = lX2=Xx + l l<yi<xi X3=X2 + 1 

xi<y3<k-2 
fc-2 fc-1 

= -T^E E E C(fc-2;[yi,xi,y3]K^^2^(^3^-4) 
1 — — — — 

Xl = lX2=2:i + l l<l/i<Xi 

= -T^E E C{k-2;[y,,x^,ysM^ ^2^(4'-^') 



xi — 1 l<yi<xi X2—X1 + I 

xi<y3<k — 2 



and since 



zj'l-^r - ^3) = fif^f^^^ ^ (22Z3)''' ^fe^r^^ - 4 



3!2=a;i + l 



Z223 ^ 1-^2 



we get 



9 fc-2 



F{k;[zi,Z2,Z3]) = -V V C(fc - 2; [yi,xi,2/3])(ziZ2Z3)' 

(l-Z3)(l-0223) ^^^^ ^^^^^^^ 
2:1 <y3 <fc'— 2 



^fc^fc+1 '^■"2 
^2 ^^3 



^ 5] C(fc-2;[2/i,xi,y3]K' 



(l-Z3)(l-^2^3),^^, l<,i<.i 

2;i<y3<fc— 2 

'7rT54r~7TE E c^(fc-2;[yi,xi,y3])(^i^2r 

^ '■^'^^ ^' 2:1 = 1 l<;^i<2:i 

2:i<y3<fc'— 2 



^fc.fe+1 *iz2 
^2 -^3 



(1 - Z3)(l - Z2) 

^ -^'^ 2:1 = 1 l<ai<2i 

xi<y3<k—2 

2 

g|2 F{k - 2; [1, ziZ2^3, 1]) 

(1 - Z3)(l - Z2Z3) 
^fc^fc+l 

(1 - 4(1 -.«) ^"'-^^ 

a^lKi'-^.) "'"'-^^"'^""'" 

^fc^fc-+i 

-F(fc-2;[l,zi,l]). 



(l-Z3)(l-22) 
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We can then use the functional equation to compute C(fc) — F{k; [1, 1, 1]) for whatever k we 
need and then find a{n) for the desired n by Equation ([S]). 

3.2. General Algorithm. In general, if we can find a functional equation for F{k,p; [zi, . . . , Zm]) 
that relates it to cluster generating functions with lower order 'fc', we can use it to compute a{n, t) 
using Equation ([S]) . One can see that most of what was done in the above example can be extended 
to any pattern (or pattern set by finding a functional equation for each pattern individually) . The 
outline of the general procedure is as follows: 

First, find the cluster recurrence for the initial summand C{k,p] [xi, . . . , Xm]) (as in Section [^?H) 
and substitute this into the summation in Equation (ITSl) . Split the summation over each summand 
C{k' ,p'] [j/i, . . . , j/m']), and handle each one separately. Rewrite the summations over xi, . . . , Xm 
and apply the finite geometric series formula as needed. Finally, express the remaining summa- 
tions as cluster tail generating functions of lower order k' . 

The only part that is not immediate is whether the summations for be ordered 

properly and whether the lower and upper bounds for each summation index can be chosen prop- 
erly so that we can adequately apply the finite geometric series formula. This can in fact always 
be done, and the ordering and choice of bounds can be done as follows: 

Let , . . . , Xi- be the entries from the original last pattern p in the length k cluster that coincide 
with entries from the new last pattern p' in the length k' cluster. In other words, Xi-^, . . . , Xi- are 
the terms that occur in the y^'s of C(fc',p'; [j/i, . . . , ym']) ■ Let Xi-^^ , . . . , Xi^ be the terms that were 
"chopped off" from the length k cluster. Also, assume that ^ < . . . < Xi. and Xi.^-^ < . . . < Xi^ . 
Note that in the example in Section l3.ll Xi-^ = xi (not "chopped") while Xir^ — X2 ("chopped") 
and Xi, = X3 ( "chopped" ) . 

Order of summations: 

The summations will be ordered (from outermost to innermost) as x^^ to Xij followed by Xij_^-^ to 
Xi^ . Thus, the outermost summation is indexed by x^^ , the next summation inward is indexed by 
Xi2 , and so on. This places the summations over Xi-^-^ , . . . , Xi^ to be on the "inside" so that they 
can be moved inward to apply the finite geometric series formula. 

Lower/Upper bounds for x^^, . . . ,Xi^, : 

For each I with j + 1 < / < m, let 6; = fc if > ij] otherwise, let hi = min({ii, . . . , . . . , «;}), 

and let q be the index of i (so ic, — hi). 

For x^^ , . . . , Xi^ . 

Xi^ = ii to k — m -\- ii 

Xi2 — Xi-^ + i2 — ii to k — m + i2 

Xij = Xi. -^ + ij — ij-i to k — m + ij 

For Xi ... : 

Lower bound is 1 if ij+i — 1, and Xi.^-^^i + 1 otherwise. Upper bound is fc — to -f ij+i if = fc, 
and — Cj+i -I- ij+i otherwise. 

For Xi, with I > j + 1: 

Lower bound is Xi, _i + 1. Upperbound is — m -|- «; if 6; = fc, and hi — ci + ii otherwise. 

One can see that the indices Xi^, . . . ,Xi. range over all necessary values and can also verify that 
Xij^j , . ■ • , Xi^ will cover all necessary values as well. Additionally, for each r, the lower and upper 
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bounds for Xi^ never depends on any Xi^ where s > r. If we applied the above approach to the 
example in Section |3.1[ we would get Xi-^ = xi going from 1 to fc — 2, Xi^ — X2 going from xi + 1 
to fc — 1, and Xi^ = going from a;2 + 1 to k. 

3.3. Additional Results. We get a couple more immediate byproducts from the algorithm in 
Section [3.21 First, the method provided for finding a functional equation always works, so we get 
the following: 

Theorem 3. Let B he a pattern set and p G B. Then, there always exists a Junctional equation 
for F{k,p; [zi, . . . , z\p\]) of the form: 

F(fc,p;[zi,...,Z|p|]) = (t-l) ^ ^ i?,.F(fe„p';[Mi,...,M|;,|]) 

p'eB iei(p') 

where I{p') is a finite index set for each p' G B, I{p') and I{p") are disjoint if p' ^ p" , each M* 
is a specific monomial in zi, . . . , z\p\, each Ri is a specific rational expression in zi, . . . , z\p\, and 
hi < k for each i. 

Additionally, we get an immediate corollary of Theorem [1] 

Corollary 1. Let B be a set of patterns we would like to avoid. Without loss of generality, assume 
that B contains no redundancies. Then by setting weight(p) = t — 1 for each p E B, we get: 



l-z~ E E F{k,p;[l,...,l])j^- 

k>l peB 

Given that we can find a functional equation given any pattern set B, in a sense, we have an 
expression for the exponential generating function P{z,t) for any pattern set. 

3.4. Maple Implementation. The algorithm from Section has also been implemented in the 
Maple package CAV. Using that algorithm, you can find the sequence a(l), . . . ,a{n) avoiding a 
set of patterns B by calling the procedure CAVT(B,n), where the patterns in B are represented 
as lists. For example, for n = 10 and B = {123,321}, try CAVT({ [1 ,2 ,3] , [3 ,2 , 1] } , 10) ; . To 
keep track of the occurrences of patterns from B, use the procedure CAVTt(B,n,t). For example, 
try CAVTt({ [1,2,3] , [3,2, 1] }, 10,t) ; . To generate the cluster tail functional equation only (en- 
coded again in a data structure that we call a scheme), use the procedure MakeTailFE(B,k,z,t) . 
For example, try MakeTailFE({ [1 ,3 ,2] } ,k,z,t) ; . 

Computationally, the algorithm in Section 13.21 is much more efficient than the one in Section 
12.11 so the CAVT procedure is much faster than the CAV procedure. In general, CAVT should be used 
instead of CAV for computing a(n) values and, similarly, CAVTt should be used instead of CAVt for 
a{n, t). 

4. Asymptotic Approximations Using CAV 

Let B = {p} be a set containing a single pattern. In [TT], Warlimont gave a conjecture on the 
asymptotics of a{n): 

(17) a{n) ~ 7 • • n! 



where 7 and p are constants depending only on the single pattern p. Some initial asymptotic 
results for a{n) were proven by Elizalde in |4j. Recently, Ehrenborg, Kitaev, and Perry prove this 
conjecture in [2]. With this result established, we can compute approximate values of 7 and p for 
various single patterns. 
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Elizalde and Noy gave some approximations of 7 and p for length 3 and a few length 4 pat- 
terns in [S]. Aldred, Atkinson, and McCaughan also gave approximations for the p values of 
the single length 4 patterns. Using the Maple package CAV, we can empirically verify these ap- 
proximations and also quickly produce many new approximations. For example, the procedure 
AsymApprox(p,N,d) will give approximate values (up to d decimal digits) for 7 and p for the 
pattern p by computing a{N — 2), a{N — 1), and a{N) and computing their ratios. For example, 
try AsymApproxC [1 ,2,4,3] ,50,20). 

To approximate 7 and p values (up to d decimal digits) for all length n patterns and then rank 
them by the size of p, use AsymApproxRajik(n,N,d). For example, AsymApproxRank(4,30,10) 
gives us the approximations for the 7 and p values for length 4 patterns: 



Pattern 


7 


P 


12 3 4 


1.1176930011 


0.9630055289 


2 4 13 


1.1375931232 


0.9577180134 


2 14 3 


1.1465405299 


0.9561742431 


13 2 4 


1.1510444988 


0.9558503134 


14 2 3 


1.1567436851 


0.9548260509 


1 3 4 2 ~ 1 4 3 2 


1.1561985648 


0.9546118344 


12 4 3 


1.1696577874 


0.9528914233 



Table 1. Approximate asmptotics for length 4 patterns 



Similarly, AsymApproxRajik(5,25,20) would give us the approximations for the 7 and p values for 
length 5 patterns. The output can be found on the paper's website. 

5. Further Work 

In this paper, we outlined the key procedures in the CAV Maple package. The cluster tail gener- 
ating function was defined, and a constructive approach was demonstrated in finding a functional 
equation for it. Using this functional equation, we were able to more quickly count permutations 
avoiding a prescribed set of patterns. In addition, by applying Theorem[2l we were able to totally 
classify c-Wilf-Equivalences in single patterns of length 3, 4, and 5 rigorously, while nearly classi- 
fying single patterns of length 6. We were also able to classify c-Wilf-Equivalences in a few cases 
of multiple pattern sets. Finally, we were able to use the faster algorithm to compute approximate 
values for asymptotic constants. 

Despite this, there is a lot of room for improvement algorithmically and quite a few new open 
problems/conjectures arise. Some of the conjectures are listed below. 

Elizalde and Noy provided the following conjecture in [S]: 

Conjecture 1. For a fixed pattern length k, the increasing pattern a ~ 12 ... k is the "maximal" 
pattern, in the sense that ao-(n) > ap{n) for all p £ Sk and all n. 

Based off of experimentation, we also have the following analogous conjectures: 

Conjecture 2. For a fixed pattern length k, the pattern a = 12 . . . (/c— 2)(fc)(fc— 1) is the "minimal" 
pattern, in the sense that ap{n) > a^in) for all p G Sk and all n. 

Conjecture 3. For a fixed pattern length k, the pattern set B — {12 . . . fc, 23 . . . kl} is the "maxi- 
mal" pattern set among sets of 2 patterns, in the sense that asin) > as'in) for all i?' G ( 2*) and 
all n. 
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Conjecture 4. For a fixed pattern length k, the pattern set B = {12 ... (k — 2){k){k — l), 12 ... {k — 

3){k — l){k){k — 2)} is the "minimal" pattern set among sets of 2 patterns, in the sense that 



Conjecture 5. For a fixed pattern length k, the pattern set B = {12 . . .k,23 . . . kl, kl2 . . . (fc — 1)} 
is the "maximal" pattern set among sets of 3 patterns, in the sense that as(n) > aB'{n) for all 
B' e (^3") and all n. 

In addition, based off of empirical evidence for single pattern avoidance up to length 6 patterns, 
we believe the following: 

Conjecture 6. For any two patterns p\ and p2 of the same length, either pi and p2 are strongly 
c- Wilf-Equivalent or they are not c- Wilf-Equivalent at all. 

This certainly holds for single patterns of length 3, 4, and 5. This also appears to hold for single 
length 6 patterns. 

Finally, the author would like to thank Doron Zeilberger for his suggestions, comments, and en- 
couragement towards the work in this paper. 
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